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An Analysis of the Effect of the Discontinuity in a 

Bifurcated Circular Guide Upon Plane 

Longitudinal Waves' 

Louis L. Bailiix* 

The problem of theoretically accounting for the scattering of sound by a semi-infinite 
circular tube with a small diameter inserted axially, as a measuring probe, into a larger tube 
of infinite length, is solved in a rigorous and explicit manner when the incident plane waves 
are restricted to the lowest propagating mode. 

The investigation follows the methods outlined by J. Schwinger and consists in solving 
Helmholtz's equation subject to boundary conditions on the acoustic velocity potential at 
rigid walls. With the aid of Green's theorem and the Green's function for a point source in 
the large guide, an integral representation of the desired velocity potential is obtained. The 
imposition of the boundary conditions leads directly to a homogeneous integral equation of 
the Wiener-Hopf type, which is solved by transform techniques. The results of this solution 
are shown to be related to a distance d, which is the length the inner guide must be extended 
to account for the disturbance occurring in the immediate vicinity of the discontinuity 
caused by the excitation of higher mode fields. 

Another satisfactory method of representing the effect of the discontinuity is by the use 
of equivalent circuits, by means of which a discontinuity in a wave guide is replaced by a 
lumped parameter network in a set of transmission lines. The equivalent circuit for this 
network is derived, and the one essential circuit parameter is related to d and to the explicit 
solution of the integral equation. 

Numerical results for d are given as a function of several different guide ratios and 
driving frequencies of the incident waves. The results are in considerable disagreement with 
approximation methods previously reported and indicate that these methods are unsatis- 
factory in the present problem. 

I. Introduction 

The pi'oblom of theoretically arcoiinting for the scattering of sound by a circular tube 
with a small diameter inserted axially into a larger tube as a measuring probe (see fig. 1) is 
merely a specialization of a general class of problems dealing with the effects of obstacles on the 
propagating modes of both acoustic and electromagnetic waves in guides. It has been shown 
by J. Schwinger [1] ^, in a sei'ies of theoretical studies, that the special class of boundary-value 



_._.i^._.- 



T^ 



2 = 

Figure 1. 

problems, where waves are incident upon a number of semi-infinite parallel metallic structures 
(structures with parallel axes) of zero thickness and perfect conductivity, can be formulated 
mathematically as an inhomogeneous integral equation of the Wiener-Hopf type [2]. Because 
of the discontinuity, the imposition of the particular boundary conditions that the field com- 
ponents must satisfy leads directly to one or more intregral equations of the type 



^W=J^ K{x-y)J{y)dy, x>0. 
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where f(x) is unknown while K(x) and g{x) are know^n functions. In the electromagnetic 
problems /(x) is the surface current density on the obstacle while, in the acoustic case, it repre- 
sents the discontinuity in the field variable on opposite sides of the obstacle. Generally g(x) 
is determined by the boundary conditions and the incident fields; it is zero in the present 
acoustic problem. The components of the field vectors at all points in space can be readily 
related to the Fourier transform oiJ(x), which is obtained when the above equation is solved 
by the transform techniques of the Wiener-Hopf method. 

Schwinger [1] also points out that if incident dominant mode waves fall upon an obstacle, 
the effect of the obstacle is to generate a far field of scattered or reflected lowest mode waves 
and a local field that is a result of the excitation of the attenuated higher modes. Although a 
complete description would include both fields in detail, as a practical guide problem, onl}^ 
the far field effect on the propagating modes need be considered. Another satisfactory method 
of representing this effect is by the use of equivalent circuits, by means of which a discontinuity 
in a wave guide is replaced by a lumped parameter network in a set of transmission lines. 

It will be the purpose of this paper to apply the above methods to an infinite circular 
guide with a concentric semi-infinite circular bifurcating cylinder and to study the effect of 
such a discontinuity upon plane longitudinal waves by replacing it by an equivalent circuit. 
By a consideration of the equations that govern them, the unknown parameter in this network 
will be shown to reduce to one single essential parameter. This parameter will then be related 
to a distance d, which is the length the inner guide must be extended to account for the dis- 
turbance occurring in the immediate vicinity of the discontinuity caused by the excitation of 
higher mode fields. By this extension, we are able to introduce a new reference plane, the 
^Wirtual plane of bifurcation.'' Here the one essential parameter vanishes, and the equivalent 
circuit reduces to three transmission lines in series as shown in figure 6. Thus all measure- 
ments made at a distance will indicate that the bifurcation occurs in the virtual plane that is d 
units from the mouth of the inner guide, provided we assume that the discontinuity produces 
no disturbance. 

The solution to the problem will be obtained by solving Helmholtz's equation subject to 
boundary conditions on the acoustic velocity potential at rigid walls. This is to be done with 
the aid of Green's theorem and the Green's function for a point source in the large guide. 
From the integrals in Green's theorem, we obtain an integral representation of the desired 
velocity potential. This integral we shall solve in its Fourier integral form with the aid of the 
homogeneous integral equation that is obtained when boundary conditions on the inner guide 
are applied to the integral representation. The homogeneous equation is of the above-men- 
tioned Wiener-Hopf type; since g(x) = Oy a solution by Fourier transforms is possible and 
yields an explicit expression for the transform of the unknown function. As the transform of 
the unknown appears in the Fourier integral form of the integral representation, the solution of 
the Wiener-Hopf equation permits a rigorous result for the velocity potential anywhere in all 
three regions (large guide, small guide, and coaxial guide). 

It is also noteworthy that an alternative integral equation formulation [3] of the problem 
can be given. This originates with the division of the space interior to the large guide into 
the two regions r>a, r<a, where a is the radius of the smaller guide. Green's theorem is then 
applied to each region with its appropriate Green's function, and each velocity potential is 
then expressed for z<CO in terms of its radial derivative on the surface r=a; here 3 is a distance 
along the guide with z=0 at the mouth of the inner guide. The requirement of continuity 
for the velocity potential on crossing this surface provides an inhomogeneous integral equation 
of the Wiener-Hopf type for the determination of the common radial derivative. 

Thus if we place the inner guide at 0>O, then, for 2<0 and large, the integral representa- 
tion reduces to linear expressions involving transmission line parameters and, for 3>0 and large, 
it is given in terms of expressions from the Wiener-Hopf solution as well as a linear combination 
of transmission line parameters appropriate to the region considered. When the solutions to 
the right and left are then transformed back to the reference plane at the mouth of the inner 
guide, we find that the discontinuity can be represented as a six terminal network and that the 
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one essential parameter and the distance d can be also expressed in terms of these Wiener-Hopf 
results. 

To simplify the problem to a considerable extent we are going to restrict the driving fre- 
quency to permit only one mode plane wave propagation in all three regions. In order to insure 
this type of propagation in the different regions, we must carefully examine the eigenvalues 
that arise in the solution of Helmholtz's equation in these regions. For the present problem 
where the normal derivative of the velocity potential vanishes at rigid walls, such a study 
will indicate that if we wish the propagating mode to be the lowest ^'axially symmetric'^ mode, 
we must have 0<^k<^Xi/b where Xi is the smallest zero of the Bessel function Ji(x). Since 
7r<^Xi = 3.832, . . ., we may therefore select w/b as a working upper limit to guarantee only 
principal-mode propagation in all regions. 

11. Description of Physical Quantities 

In studying the effect of the discontinuity upon a plane acoustic wave in an infinite circular 
guide, radius b, which contains a concentric semi-infinite coaxial guide of negligible thickness, 
radius a, we will assume that all the walls are perfectly rigid. Accordingly, the radial com- 
ponent of the particle velocity vanishes on the rigid walls. Thus, the scalar velocity potential 
* satisfies the wave equation 

''^=c'V% . (1) 

where c is the velocity of sound propagation in free space, and the partial derivatives of $ with 
respect to the coordinates are the respective components of the velocity. A solution under the 
assumption of harmonic time dependence and complete angular symmetr}'' about the 2:-axis 
can be given as 

Hr,z,t) = ct>{r,2)e-'^'', (2) 

where 

and X=free space wavelength. Equation (1) then becomes the scalar Helmholtz equation 

(v^+^^)*=0+7§^+||+^-V=o, (3) 



subject to the boundary condition 

(r=b for all z 

(^> 

[r=a for 2>0. 



^=0 at f 



Once is known, we can calculate the other important acoustic quantities from the classic 
relationships between the velocity potential, pressure, density of gas, and condensation. 

We now formulate the equation which expresses the field variable in terms of the dis- 
continuitj^ of across the surface of the inner guide (surface current density in electromagnetic 
case) . To do this, we are going to make use of Green's theorem to construct the Green's func- 
tion for the structure at hand out of a more readily obtainable one, which in our problem is 
the one for the empty guide of radius b. Thus we are considering the inner guide as an obstacle 
embedded in the space of the other. The appropriate Green's function which will aid appreci- 
ably in the solution of the differential equation in (j) satisfies a similar but inhomogeneous 
equation with the simpler boundary condition 

-^=0 at r=b for all z. (5) 
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Therefore, Green's function for the velocity potential, with no angular dependence, satisfies 
the equation 



Here b{x—x^) is the Dirac delta function and is zero everywhere save at x—x\ where it becomes 
infinite in such a fashion as to make 



p b{x-x')dx'=\. 



(7) 
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Figure 2. 







In applying Green's theorem,^we take the region under consideration as that which is 
enclosed by the dotted Une in figure 2 and write 



J^^ [GiV^+k^.i,-4>iV'^+k^G]dV'=l{6 ^-* g) dS', 



(8) 



where n'= outward normal to the bounding siu-face. Although dV^ represents an element of 
volume enclosed by the dotted line and dS' an element of surface area, the problem can be 
reduced to two dimensions by immediately integrating with respect to the angle, since there is 
complete circular symmetry. This integration merely brings in a factor of 27r on both sides of 
the above and removes the angle dependence from all subsequent equations. G here is an 
appropriate Green's function that satisfies (6) and (5) . 

To evaluate (8), we impose boundary conditions (4) and (5) together with (3) and (6). 
Since the volume integral reduces to <^(r,s), Green's identity yields 



♦('•^'"K" »-*»)"«'. 



(9) 



where S^ is determined by the six regions in figure 2 and n' measures distance along the out- 
ward normal to the surface enclosed by dotted line. By writing 6~ir, r\ z—z') as the asymptotic 
form of Green's function for z' large and negative and G'^(r, r', z—z') as the asymptotic form for 
z' large and positive, (9) becomes the sum of the six integrals. 
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(minus in front since n'=—r'), 

(minus in front since ri' = — s'). 

Since it will be shown later that all integrals taken along transverse surfaces are independent 
of L or U as the case may be, we may now let L and U -^<=o and write (10) as 

*<-)-r<*^-«(a..«"''^L"i[j::(''*i'-*i^),. 
^-J:(-l?--m.'--j::(-.^^4!").....-] <■•> 

The differ(\nce of the velocity potential on the inner and outer surfaces of the small guide which 
appears in the first integral of (11) is a consequence of the oppositely directed normals at 
these surfaces. 

III. Derivation of the Appropriate Green's Function 

To obtain the solution G(ry ^z) of (6) in the total space interior to the larger guide, we 
shall apply standard Fourier transform techniques. Thus, we transform z out of (6) by de- 
fining the Fourier transform in z as 

2^(r,r',f)= f"^ e~''^G{r/,z)dz, (12) 

»/ — CO 

where we have written z instead of z—z'^ since this binomial appc^ars explicitly in the solution, 
^ now satisfies the equation 

a"+;^+^')*<-''.»==^, (.3, 

where 7^=F— f^;2^is subject to the boundary conditions: (1) that it be finite at the origin 
and (2) that its normal derivative vanish at r=h. Thus, we must seek a solution of BessePs 
equation of order zero, since there is angular symmetry, which satisfi(^s thesc^ boundary con- 
ditions and is continuous at r=r'. Following Courant-Hilbert [4], the solution of (13) can 
be given as 



5^^(^ r', t)=\ ^"^j \Jlir)Nhh)-Nlir)J,{'ih)\ {r>r') 



(14) 



By the Mellin inversion, these then yield G in the form 

^^iX-T '^'''"^^ ^J^ [ Jo(X/'>)iVi(7^')-iVo(7r>) J,(t6)] (1 5) 
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where the integration contour is a straight hne in the region of regularity of the Green's func- 
tion transform and r>, r < symboHzes the larger and smaller of the coordinates r, r', respectively. 
To determine the region where the integral in (15) converges, we must first determine 
the form of G and then apply it to (12) which defines ^(f). From our fundamental assumption 
of only one mode propagation, we know that 



G^ 



71 = 1 



(16) 



It can be easily shown that for a 6* of this form, the integral in (12) has no region of regularity 
unless k, the propagation constant, has an arbitrarily small positive imaginary part (which is 
eventually set equal to zero). This corresponds to a small attenuation of the sound waves 
traveling in the guide. Thus, if we let k=ki +ie and introduce the coordinates in the f-plane 
as f=S+ ii], (12); converges in the strip |i7|<e (Im k). Consequent!}', the inversion given 
in (15) can be performed, if a is restricted to this same strip. 
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Figure 3. 



t^- PLANE 

The integral in (15) may now be evaluated by computing the residues inside the semi- 
circular contours pictured in figure 3. The poles of (15) which are given by the zeros of Ji{'yb) 
are all complex for e>0. Thus, if we let Jj(7„6)=e/i(a:J = for 7i=0, 1,2,3, . . . and ^n= 
k^—yn, compute the residues inside Fj and F2, and then let e— >0, we find that the solution of (6) 
becomes 



G= 



ieik\z-z'\ 



kb' 



■s 



Knb^JliXn) 



Here we have insured only one mode propagation by restricting k to be less than Xifb and 
defining the attenuation constants by ¥.l={XnlbY—¥ where K\={xilby—k'^ is the smallest. 
The Green's function which is appropriate to this problem must give a physical descrip- 
tion of waves incident upon and reflected from a discontinuity in a guide. Since (17) repre- 
sents a wave traveling outward from a source, the appropriate Green's function requires a 
superposition of (17) and its complex conjugate which is the solution of (6) representing a wave 
traveling inward to a sink. Thus if we take one half of the sum of (17) and its conjugate and 
simplify the results, we obtain 

a,.,. .,,_ sinA.|.-.1 - '^°(--i)'^°(--0 --""'---'' 

Cr(/,r ,^ — <^ ) — z;:7;-2 rZ-! Z^~1J1^7Z~\ ' (1^) 



kb' 



which has the desired physical properties. 



K^b^Jl{Xn) 
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IV. Evaluation of Integrals 

Once the appropriate Green's function has been (leterniined, we can proceed to evaluate 
the integrals in (11). This is done in paragraphs (a) and (b), since the character of the results 
and the evaluation techniques used differ widely. Paragraph (a) is concerned with obtaining 
the far iield contribution in terms of the transmission line parameters and paragi-aph (b) with 
the setting up of the Wiener-Hopf integral equation which is a result of the first integral in (U). 

(a) In proceeding, let us temporarily ignore the first integral of (11) and examine the last 
three. From the explicit form of the Greenes function, (18), the asymptotic forms of G can be 
given as 

G^{r,r',z-z') = G-{r',r,z-z')=- ^'''^l\-''y (19) 

Thus, when we examine (11) for L and U sufficiently large so that the attenuated terms may be 
considered negligible, only the forms in (19) need be used. It can be shown [5] that in terms of 
the propagation of acoustic waves in transmission tubes, the velocity potential and j)ressure 
obey analogous equations to those for currents and voltages in a transmission line. Thus, for 
z large, </> behaves as the current in a transmission line which we will denote as I{2). Then, in 
terms of the transmission line currents and voltages measured at 2 = 0, we write 

<t>(j% z)^I{z) = [I cos kz — iVY^m kz], (20) 

where 

z pc 
is the characteristic admittance of the line and p is the density of the air or the gas in the tubes. 
Using (20) for each of the three regions under consideration, we note that as s-^— co 

(t>(r,z)'^I,(z) = [Ii cos kz+iViYi sin kz], (21) 

while for z^<x> and for r<ia 

<t>(r,z)^U2) = [I, cos kz-iV,Y, sin kz], (22) 

and for 2->co and for 6>r>a 

(I>{r,z)^l2{z) = [l2 cos kz-iV^Yo sin kz]. (23) 

Thus, the last integral of (11) becomes 

«ii[-/:(«-i?---')-*m=-."*']-t-*' 

since 

d</) J. r • 7 /, T/T^ 1 f^ n- ^\Vik{z — z') bG' 1 J. ,, 

^— ,=/:(— ii sm kz +1 Vii i cos kz ), G = yu ^ ^r~r=-r^ cos k{z — z), 

oz ' Kb" Oz b~ 

Likewise, 
and 



since 



^=k {- 1, s^m kz-iV.Yi cos kz'), (?+ = ^, — ', -^= ^ '- 
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Thus (11) becomes 

<^(r,.) = J>.')(gX^/.'+^[/.(3) + (^)/.(^)+(|;)/3(3)} (27) 

where 

H{z) = a [0+ {r,z)r-^ai.o-^~{r,z)r=a-o\' (28) 

Since (28) is independent of L and U we, therefore, note that all the integrals taken along 
transverse surfaces yield propagating terms which depend only on z. Equation (28) represents 
the velocity potential at all points in the large guide in terms of the potentials of the three 
far fields and the potential of the unknown scattered field which is given by the integral. 

(b) Since H(z) is as yet an unknown function, (27), the equation we referred to in our 
introduction as the integral representation of the desired velocity potential, is an inhomo- 
geneous integral equation. When we impose the acoustic boundary condition 






0; for3>0; 



this equation becomes readily soluble by transform methods since it is homogeneous and of 
the convolution type. 
Thus, if we call 

and since the propagating terms drop out in the differentiation, we have from (27) 

0=r K{z-z')H{z')dz', 2>0 

V{z)= r K{z-z')H{z')dz', 2<0 

where the range of integration has been extended to — <» , since, from (28) 

(0, s<0 



(29) 



H{z)-- 

{H(z), 2>0. 

Equation (29) is a homogeneous integral equation that resembles Wiener-Hopf type, since 
the dependence of G and therefore the kernel, 



^(^-^')-KP^'->'^-^'^1Lj 



upon z is of the form (z—z^). 

V. Derivation of the Lumped Parameter Circuit and the One Essential Parameter 

From equation (27), we can now derive the circuit relation of a lumped parameter net- 
work with which the discontinuity can be replaced, if we concern ourselves with fields far 
enough away from the mouth of the inner guide so that all attenuated modes have essentially 
vanished. To investigate this far field result on the left as z-^— oo ^ we let (l>(r,z) be given by 
the transmission line expression (21) and define the characteristic admittances as inversely 
proportional to the area: 

11 111 

(30) 
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Thus, as we shall see later, the integral in (27) vanishes for z<iO and large and we are left 
with 



which becomes 



ZJ,{z)+ZJ,{z)-ZJ,{z) = 0. 



(31) 



Now, if wc rewrite the above in terms of the /'s and V's at 2=0, we have 

Z2(/2 COS kz—iV^Yi sin kz)-\-Zi{Ii cos kz—iVzYz sin kz) — Zi(Ii cos kz-\-iViYi sin kz) = 0, 



and since the sine and cosine are Hnearly independent functions, (31) becomes 

ZJ,+ZJ,-ZJ, = 

F,+ F3+F3=0. 



and 

These two, together with 



(32) 
(33) 

(34) 



a result of (30), are the defining equations for a six terminal lumped parameter circuit. Thus 
vve have three relationships between the circuit parameters which were deduced by examining 
the effect on the velocity potential as 2— >— 00 and which represent the far field results to the 
left. We will presently determine the far field results to the right and then by expressing 
them in terms of the current and voltages at the zero reference plane, we will find that the 
extrapolated results coming from the left and right can be represented by same six terminal 
impedance network. 

In general, a six terminal network contains ninc^ impedance parameters representing the 
ratios of voltages to currents for all possible combinations of the three independent voltages 
and currents. However, if we make use of the reciprocity theorem (see Schwinger [1] or 
Everett [6]) of electrical circuit theory the number of impedance parameters is reduced to six. 
From the above three equations, we further reduce this number to three, which we will call the 
susceptances, Bi, B2, and B^. Also, from (34), we notice that they must be proportional to 
Yi, Y2, and Ys, respectively. If this circuit is to satisfy the above equations, it must have the 
following delta form (fig. 4) with nondissipative pure reactive elements and with the power 
flow arbitrarily chosen to flow into the junction from both sides. 
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Since the circuit parameters are independent of the impressed currents and voltages, we 
can now obtain independent relationships between the i^'s as follows : 

1. Open circuit lines 2 and 3. Thus l2=I^=0 and from (32), /i=0. By looking into 
line 1, we see Bi in parallel with the series combination of B2 and B^ and we have 



KBi- 



iB2^iB,) 



= or 
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\_Bi B2 SgJ 



^0. 



(35) 



2. If we terminate lines 2 and 3 by their characteristic impedance, it can be easily shown 
that the impedance seen from line 1 is the characteristic impedance 



Y^=iB,+ 



{iB,+ Y,){iB,+ Y,) 
{iBri-iB,+ Y2+ F3) 



(36) 



With the independent derivation of these last two equations, we see that the number of 
independent parameters is now one. Thus from the proportionality required by (34), we let 



B,=ba\, B,=-.h2Y„ and B,= h,Y,. 



Then from (36), wo obtain 



and (35) becomes 



Y,{\-ih,) Y^{ib,+ \yY,{ih^^-\) 



(37) 



(38) 



By solving (37), (38), and (34) simultaneously for the three unknown parameters, we find that 

bi = -b2=-b3=B (39) 



and figure 4 reduces to figi 
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Thus we have shown that with respect to the far field results on the left, the discontinuity 
can be replaced by the above network. To determine the effect on the right as s->+ oo, we 
must first evaluate the integral in (27). A straight forward circuit analysis of figure 5 shows 
that B is related to the circuit parameters b}^ 



B= 



iih-h) 



F1F1+F2F2 



(40) 



We have now shown that the circuit in figure 5 can be used to represent the discontinuity, 
when the parameters are measured from the ^=^0 reference plane which is at the mouth of the 
inner guide. Since, there is only one essential circuit parameter which has a known variation 
with driving frequency and the distance to the reference, it is possible to find a reference plane 
in which the equivalent circuit susceptances vanish and, consequently, the equivalent circuit 
reduces to three transmission lines in series as shown in figure 6. 

This new reference plane is a 'Virtual plane of bifurcation", inasmuch as all measurements 
made at a distance will indicate that the bifurcation occurs in this plane, provided we neglect 
the effect of the discontinuity. Particularly, if one of the tubes is introduced in order to take 
measurements, the quantities will be measured in the virtual plane and not at the mouth of 
the tube. 

To find this virtual plane which is at 2= —c?, or a distance d ^ 'ahead" of the actual plane 
of bifurcation, we follow the method outlined by Miles [7] and let rj^^kd and transform the 
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^ 2 



^T', 






Figure 6. 



impedance seen by the i-th tube (i= 1,2,3,), Z^. (plane z=0), to ZV^ (plane z=—d). This 
transformation is effected by the use of 



2' =^~7^ ^^^ ^+i sin 7] 






± cos y\-\-i 



(^')sm. 



(41) 



hut car(* must be taken to select the ])i'oper sit>;ns (botli plus if the direction of the* power flow 
obtained by following Vi then li is in the positive z direction) and to use a negative* t] since 
the reference plane is being moved to z=--—d. Thus, if we load terminals 2 and 3 (fig. 5) with 
arbitrary athnittancc^s Yt^ and Yt^, it can be easily shown thnt liri(» 1 sees 



Thus if 



tluMi fi-om (42) 



_ 1 _. , [YT^-iBY2][YT^-iBY^\ 

^^^=z;:-'^^'^+ [Yr,+Yr,-iB{Y;+n] ^^^^ 

2^3-0 = -.^ and Fr^^O, (43) 



^n= iBiY._Y.) ^^^) 



Now, if we transform the above Z^. to Z^. until the circuit can be given as the scries com- 
bination (fig. 6), namely, 

z;,=z;,+z',3, (45) 

we find that 

i?=tan r). 

VI. Fourier Transform Solution of the Integral Equation 

The Fourier transform solution of (29) will be undertaken as follows: In part (a) of the 
present section, the equation linking the transforms of the respective functions y, //, and K 
and the common region of regularity of these transforms will be (kM-ived. The equation follow^s 
as a direct consequence of the convolution theorem. The regions of analyticity of the trans- 
forms are obtained from the asymptotic forms of the functions themselves. Part (b) will 
concern itself with the decomposition of the transform of the kernel, K, into the quotient of 
two functions which will convert the equation of the transforms into a single integral function. 
This integral function will, in part (c), be shown to be constant since it satisfies the conditions 
of the Liouville theorem. Thus from (70), we will be able to find an expression for the trans- 
form of the heretofore unknown function H{z)y and thereby, in part (d), we will obtain <^(r,s) 
from the inversion integral. 
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(a) The solution of the extended Wiener-^Hopf integral equation (29) by Fourier transform 
methods can be facilitated by the theorem of convolution which permits us to write the trans- 
form of the product of two functions as the product of the transforms. Hence, the transform 
of (29) becomes 

yx(f)^^^(^)jr(r), (47) 

where 



jr(r)= 



r-s-''-[um....j^'-Hm-)i,.. <-' 



and yf{^) and "/^(f) denote the Fourier transforms of the respective functions, H(2) and V{2) 
in the form /(f) = j{x)e~^^^dx. For the transform equation (47) to have significance, there 

J -00 

must be a common domain of analyticity of the transforms in the f -plane. To verify the 
existence and location of this domain, we examine each transform with the aid of the growth 
order of their corresponding functions V{z), H(z) and K(z). 

Previously, we showed that ^'(f) exists only if k has an imaginary part. However, >f^(f) 
involves differentiation with respect to r and r' which eliminates the propagating term and, 
consequently its region of regularity is |i]| <Ki. 

The quantity which we have called H{2) , given explicitly by 



[(t>{r,z)r=a-hO—(t>(r,z)r=a-o\(^ 

has the same asymptotic form, for z-^co^ as the velocity potential on the inner and outer sur- 
faces of the smaller guide and, therefore, can be given as H{z)-^Ae'^^-\-Be~'^^. From this, we 
infer that Jf'(f) has only simple poles on the real axis at f = ± A:, and is regular in the lower half 
plane ?7<C0, since 



is bounded in this region. 



\ynn\<£mz)y'dz. 



(49) 



If we study the defining equation for V(z), (29), by inserting the explicit expression for 
the asymptotic behavior of the kernel, namely K{z— z^) ^e~^^^^~^'^ , we see that the real part 
of V{z)^e^^^ as 0— >— 00^ since H{z) is the entire function (traveling wave) given above. 



Thus 



and exists for ry^ — Ki. 



f"(r)~£^ 



g-ii-j+/<r,2^2, 



-L 



e-i^^+(v+K,)zdz 



(50) 



>^CS) 



;h:(sk 



L L 



^{t) 






>-COMMON REGION 



■S - PLANE 



Figure 7. 
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Collecting results, we find (fig. 7) that tlie strip —Ki<Cr)<CO is a common region of regu- 
larity for the transform equation (47). It will now be convenient to designate the region 
77> — Ki and rj<iO as the upper and lower half planes respectively. 

(b) In order to solve the transform equation (47), we now seek to represent .vr(f) as a 
quotient of functions K_(^) and K+(^) such that a simple transposition will make each side of 
the equation regular and not zero in separate half planes with a common strip. By writing (47) 
as 

7^(^)7?+ = jeU)K_ = integral function, (5 1 ) 

we will see that each side of the equation is the analytic continuation of the other and there- 
fore will represent the same function. Together they will then define an analytic function 
throughout the finite f-plane. If, in addition, this analytic integral function turns out to 
be bounded in the whole plane, then, as a consequence of the Liouville theorem [8], it must 
be equal to a constant and the transforms ,^(f) and ^(f) are simply related to i^4- and K^. 
Although there are general methods for the decomposition of the kernel which are both 
rigorous and forceful and have been successfully [3, 9] used, in our present discussion however, 
a simpler and more direct approach will suffice as it has in other problems [10]. Thus we 
consider 

and write each factor as an infinite product of the form {(i^—¥) which we will simply split 
into terms which have simple zeros in one half plane or the other. 

Thus to convert (52) to infinite products, we start by letting 

j{y)=J,{ya)N,{yh)~J,{yh)N,{ya), (53) 

From Whittaker and Watson [11], this may also be written as 

m=cn (\--l'-\ (54) 

where yn = i_ +Q(~) ^^'^ roots of f{yn) = 0, provided that 

(1) f(y) is an even function 

(2)/(0) = C' 

(3)/'(0) = 0. 

A careful analysis of these restrictions will show that they are satisfied, if c=^{b^—a^)/{irab). 
Thus, we can now write 



where e ^"^ are inserted as absolute convergence factors, without which the infinite products 
are only conditionally convergent [12]. Similarly, if we expand 
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where Jn{Xn)^0, then from (56) 



a^'k^-^' 



and 



mr 



Hence, we may have 



(57) 



Jf(f) = 



i^ib-a) 



n!V^'s^l« "Sjv. 






jnb-a) 



[5!V'-f-¥l-v!V-f-^!-l 



(58) 



(c) Now, if we define, as shown in (58) , 



KM) 

K4n' 



(59) 



where K^) is regular and not zero m the lower half place, ?7<0, and K^) is regular and not 
zero in the upper half plane, 7?> — ei> — Ki, our transform equation (47) becomes (51) which 
satisfies the requirements of the Liouville theorem because '>^(f)jK'+(f) is regular for r?> — €i, and 
Jt{^)KM) is regular for 77<CO and the strip 0>7;> — ei>Ki preserves the continuity. 

We must now study the growth order of both members of (51) as |f|->oo. Thus for 
f {^) and J^(f ), we examine the transforms 



^•(f)=j; 



e-'''H(z)dz. 



(60) 



(61) 



The behavior of F(2)= -^ (r, z) as 2^0, can be obtained by examining 4> in the region of 

the sharp edge (zero thickness). In this region, where the values of z are very much smaller 
than a wave length, <j) satisfies Laplace's equation w^hich can be solved by static methods. 
This is done by transforming the potential field about the sharp edge to the uniform region in 
the upper half plane with the aid of Schwarz-Christoffel [13] transformation, W^=2z. From 
this transformation, W ^z^'^ and the field which corresponds to V{z) as s-^O at the sharp edge 
is \dWldz\^2z~^^^. Therefore, V{z) which is related to the component of the particle velocity 
along the c^dindrical radius at the inner pipe has an integrable singularity of the form z"'^'^ 
and 



J -00 Z 



(62) 
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Also from the above transformation, we see that H{z) which behaves as the velocity potential 
itself can be given by H{z)^z^'^ as ^-^0 and 



yn^)- ^^ e-^^\z"'')dz^ j^,^^ lf|->- for 77<0. 



(63) 



Thus, from a consideration of the physical nature of V and H, we find that Jt and / behave 
in a purely algebraic manner as |f |— >oo. Therefore, the acceptable results of the decomposi- 
tion of the kernel must behave in a similar fashion to permit the integral function to be 
bounded. 

To properly decompose the kernel, we choose 

where again 



ynir J . \nnj 



, and e 



are insertinl in each case as absolute convergence factors and the arbitrary factor e^^^^ is put 
in to remove the (exponential growth and thereby insuring the algebraic l)ehavior of the 
X's by the proper choice of x(f). This choice will be determined from the asymptotic form 
ofX_(f)[orK+(f)] as Ifl-^cx. for r7<0. 

We then find that K^, has an infinite niunber of poles on tlie negative imaginary axis for 



'<-^|-'-(J)\D—'^>--^■ ^«4> 



and zeros for 77<— Kiand is therefore regular and not yato hi the upj)er half ])lane r;>— f^i. 
By removing (63) from (58) , we have left 






^x (f) 



7i:_(^)^ ^ ^ /x ^ — — ^ : _!: ■ ,/ J — — ^ (^'^) 



which is regular and not zero in the low^er half plane, 'n<CO, and has only simple poles and 
zeros in the upper half plane r;>0. We note, however, that K^ and K^ have a common 
region which is at least as large as 6]. 

For the purpose of evaluating the asymptotic behaviors of K- and K.^ we may write 
yn^mrl(b — a) and, since Schwinger [1] proved that the infinite products in (65) must be inde- 
pendent of k as If |->oo, the square roots behave essentially as unity. Therefore, as |f|->oo^ 
we may use the Weirstrass [14] definition of the V function and the Sterling asymptotic for- 
mula [15] for r(f) to write, as a typical example, 



fi!V 



1-^+-^ e «' ~ ==.» |rH°°, ')<0. (66) 



where a: = ^f(& — oO/tt aiul (7=EuIer's constant. Similarly, 



u 



— -L '•e «'~ ' ' , (67) 

gC.g-(.+i,(y^i)i^2;r(2/ + i) 
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i^a 



where y = 

Thus as If 1^0, we have 



m 



-yj2i^be^ 



(D 



r) 



■yl2i^a(b — a) 



C"^) 



(68) 



V2if(6-a) 



which approaches (if )''^, provides that 



-ir 



(69) 



x(f)=^^ [b lnb-alna-{b-a)lnib-a)] 

IT 

to insure the algebraic behavior of both K- and K+. Likewise, we may show that 

ii:+(f)->(-^f)^ or that jr(f)->f. 

Collecting results we see that "/^ (f) K^.(^) is bounded as |f|-><» for rj^—ei and ^(T) 
K-{^) is bounded as |r|-^«> for i7<0. Thus, we have shown that both sides of (51) define an 
analytic function which is regular everywhere including a common strip of width 6i and is 
bounded at infinity. Therefore, by the statement of the Liouville theorem, (51) becomes 



T^(f)K+(r)= jeit)KM)=c. 



(70) 



If we were interested in the difference of the velocity potential on either side of the inner 
pipe, we could now obtain it by evaluating 






KM) ' 



where Br is an infinite contour within the strip of regularity |r]|<K:. (d) We are now in a 
position to evaluate the Fourier integral form of (27), which we may write as 



^Tj!,+= \_Or J,'=a, 



(71) 



and where we recall ^ exists in the strip— Ki<^r]<CKi and JffCf) exists for r;<0. Therefore 
the path of integration must be taken in the region 0>rj>— Ki. 
If we concentrate for the moment on r<^r', we have from (14) 






=^ y #?; [- J^iya)N,iyb) + mya)J,{yb)] 



Jiiyb) 



and 



^(f)=f y' ^ll^ [Mya)myb)-N^{ya)J,{yb)] ^ 
KM) 



KM) 
Therefore, if we form the ratio 



(72) 



(73) 






Joiyr) 



then 



^(f) 


yJiiya) 


dG 


Joiyr) 


KM) 


dr' 


yJha)KM) 
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(74) 



(75) 



Now the integral in (71) becomes 

which has simple poles at f=±A:(7 = 0) and f=itiK^. 

The contributions at the poles on the imaginary axis vanish as |s|-^cx3 since they contain 
g-i^„|2| Yov 2<C0, we closed the path of integration by an infinite contour in the lower half 
plane. Since the poles enclosed are all on the imaginary axis, the contribution to the integral 
vanishes as 3->— oo and the result claimed in section V is established. For s>0 and large^ 
we use a contour in the upper half plane which encloses f=±^ and (71) becomes 

since only the residues at the poles on the real axis contribute as z^>co . Likewise, if we con- 
sider ry>/, it can easily be verified that the result is the same as (77) with a different constant 
D. 

VII. Expression ''d'' in Terms of Infinite Series 

In (77), we have evaluated the integral representation (27) from its Fourier Integral 
form (71) and we can now obtain the far fields to the right. For now, as 3— >oo^ 0(r, z)-^I^{z) 
in region 3 and 12(2) in region 2. In terms of the parameters at z—0, we have similar results 
in regions 2 or I] with different constants A namely 

- % 5^ [ F3 r,+ V, Y,- V, F,]. (78) 

If we, again, equate coefficients of sin and cos, we have the far fiekls results extrapolated from 
the right as 

and 

[-iF.,3F,3+i(F3F,+ F.F,-F.F,)]=A3[^f=^-^} (80) 

Recalling that the results from the left gave the circuit equations (32) and (33), we apply 
them to the above and find that 

[/,,3-/,] = ^2,3[^^(_^p^^J (81) 

and 

-^(F.3F,3+F.F,) = -i^.3[^y+^} (82) 

By simply dividing these two equations, we see from (40), that the left side is merely the one 
essential circuit parameter 6 discussed in section V and hence with (46), w^e have 
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where the reciprocal of K+{k) is 



-( L k'a' ika\ ^ -( r k 

1 ?W^-^-^r ?w^-T 



a;^ a; 



?{V^ 



7„ Tn) _ ^g^^ 



X (A;) 



B^^ writing the bracketed terms in (84) as exponentials, for these are of the form cos d—i sin 6, 
we see that (83) becomes 

f =[>„ „<i-Oi«<»- .)Hshxs)-fe)4s[--<5^)-(^)]+il;[™ii:)-©} 

(85) 
where yn are the roots of 

JiiVn) N, {ayn) - J I {ayn) N, {yn) - (86) 

and Xn are roots oi Ji{Xn) = 0. For ease in computation, we have made the following parameter 
changes: 

1. Let a^hja. 

2. From our initial restriction on the driving wavelength, let k = rj7rlb, where 0< r)<l. 

3. Let {yn(i)=yn' 

Since we now have an analytic expression for d, it is of considerable interest to determine 
^ as a function of frequency and the radii of the tubes. To evaluate d numerically, mustwe 
find some way to get an accurate result for the above infinite series which is of the form 

± {[sin-H^n)-F,]} = i: {[sm-KXn)-{Xn) + {Xr.-Y,)}, (87) 

71=1 « = 1 

The first two terms can be written as 

f:[sm-\Xn)~Xr,]+Rn, (88) 

71 = 1 

where N is some arbitrary finite number. Thus by calculating A^ simple values and bounding 
the remainder, we will be able to evaluate these sums Writing sin~^ as an infinite series, 



Rn = ±lsin-\X„)-X,] = ±^+l(X„)^+f^{X„y+-^{X„y+ . . .] 



(89) 



we have a series which converges rapidly for XnK.1. 
If we now define Un+i as 

TT sin-^(Z^+i)-X„^.i ^ 1 

we see that the numerator is the largest term in the remainder and that U^T^l approaches 1/6 
as N grows large. Now from the series expansion, (89), it can be easily verified that 

sin-^XT^+i) — Xjv+i sin~HXjv+i) — Xjv+i , 

{x^^,r - {x^^^y ' ^^'^ 

for i—l, 2, 3, since Xiv+i>XAr+i. 
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Thus from (90) and (91) 




U^+AXy+,y> sin-' {X^+,)-X^+t 


(92) 


which means that 




Rn<U^+^^^iXr,Y. 


(93) 


Also from the series, we have, by dropping higher order terms. 





A-1-1 O 



which to^rother with (93) bounds JS„ as 

1 



O iV-f 1 A^+1 



(94) 



(95) 



Tims, we can easily determine a numerical value for R^. The remaining sums of (87) may be 
easily obtained since lim (X„)= F^. 



The infinite sums in (85) were evaluated by the computation laboratory of the Institute 
for Numerical Analysis, National Bureau of ^Standards, for r? , the frequency parameter, assuming 
the values of 1, .8, .6, .4, .2 and 0. Although 77 = has no physical meaning in the present 
acoustic problem, for completeness the mathematical lim ird/b was evaluated to show that d 



rj-^O 



is a smooth and continuous function of 97 for 17 ==0. Thus shice 



equation (85) becomes 



lim sin ^ 7^=7t? 



T=["' "+(-!-') '"<"-"]- 






^n 



+ 



-1 

aynj 



(97) 



as Tj-^0. The values of a used were 1.05, 1.2, 1.5, 2, 8, 5, and 10, whih* a^ 1.0 and 00 would 
mean that there is no inncu- guide; consequently d = 0. The computations were intended to 
be accurate to four decimal places, but with the methods outlined greater accuracy could have 
been obtained. Tlu^ eigenvalues yn (table 1) for small n were obtained from Muskat et al 
[16] with corrections given by Miller [17] and for large // were computed with the aid of a 
formula give by McMahon [18]. The results of (85) are exhibited in tables 1 and 2 and in 
figures 7 and 8. 

Table 1 











Vn 








Xn 


\ 


















1.05 


1.2 


1.5 


2.0 


3.0 


5.0 


10.0 




\ 


















1 


62. 837 


15.728 


6. 3218 


3. 1966 


1. 6356 


0. 84715 


0. 39409 


3.8317 


2 


125. 666 


31.426 


12. 586 


6. 3124 


3. 1789 


1.6111 


. 7330f) 


7.0156 


8 


188. 497 


47. 131 


18. 863 


9. 4446 


4. 7381 


2. 3853 


1.9748 


10.1735 


4 


251. 329 


62. 837 


25. 143 


12. 581 


6. 3028 


3. 1642 


1.4189 


13. 324 


5 


314. 160 


78. 544 


31. 424 


15. 720 


7. 8696 


3. 9454 


1.7643 


16.471 


6 


376. 922 


94. 251 


37. 706 


18. 8()() 


9. 4375 


4. 7279 


2.1107 


19.616 


/ 


439. 824 


109. 958 


43. 988 


21. 999 


11.007 


5.5112 


2. 4578 


22. 760 


8 




125. 666 


50. 270 


25. 141 


12. 576 


6. 2950 


2. 8052 


25. 904 


9 






56. 553 


28. 281 


14. 146 


7. 0791 


3.1529 


29. 047 


10 






62. 836 


31. 422 


15. 716 


7. 8635 


3. 5010 


32. 190 


11 








34. 563 


17. 286 


8. 6480 


3. 8492 


35. 332 
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9. 4327 


4. 1975 


38. 475 
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10.217 


4. 5459 


41.617 
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4. 8944 


44. 759 
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5. 2430 


47. 902 



333 



Table 2 
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1.05 
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.4802 
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1.2 
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1.5 
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.5740 
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.2000 
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.1752 
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1 1 1 
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Figure 9. 



VIII. Conclusion 

In (85), we have obtained a rigorous and explicit solution for the parameter d. An approx- 
imate solution of the above problem was first reported by J. W. Miles [7], who showed by 
using impedance concepts developed in earlier papers and the variational principle of Schwinger's 
that the equivalent circuit elements can be related to the analogous change of cross section. 
Thus by considering the bifurcation as a mere change of cross section, he obtained his equation 
(53) which contained only the leading term of an infinite series. A comparison of the results 
given in table 2 and the values given by Miles indicates, that for this problem, this type of 
approximation is very poor. Attempts have been made to improve this method of approxi- 
mation by considering two terms of the infinite series and only shghtly better results were 
obtained. An explanation of these discrepancies between the exact results given here and 
the approximations is to be given by Miles [19]. 



The author is deeply indebted to Professor Alfredo Banos, Jr., for the guidance that he 
has so graciously given throughout all phases of this research and to Professors Banos and 
David S. Saxon for the suggestion of this problem, as well as the benefit of their experience 
in the Radiation Laboratory of Massachusetts Institute of Technology. 
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